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Abstract 

We study the effect of scalar leptoquarks on the exclusive rare B meson decays 
K*^ K-7r+)l+l- in the full kinematically accessible physical region. We work out the con¬ 
straints on leptoquark parameter space using the measured branching ratio of Bg —)■ process 

by the CMS and LHCb collaborations. We compute the branching ratio, forward-backward asym¬ 
metry and isospin asymmetry distribution using the constrained parameter space. We also look 
into various form factor independent and CP violating observables in the scalar leptoquark model. 
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I. INTRODUCTION 


The study of rare flavour changing neutral current (FCNC) transitions of 6-flavored 
mesons decaying into dileptons provide an ideal testing ground to critically test the stan¬ 
dard model (SM) and to look for the possible existence of new physics (NP). Such pro¬ 
cesses are highly suppressed in the standard model as they proceed through amplitudes 
involving electroweak loop (penguin and box) diagrams. Of particular importance are the 
rare semileptonic decays involving b transitions, as these processes are one-loop 

suppressed in the SM, but many extensions of the SM are capable of producing measur¬ 
able effects in various observables. While most of the flavor observables are in very good 
agreement with their SM predictions there are some exceptions in semileptonic B decays. 
Recently LHCb has reported deviations from the SM expectations in R —)■ angular 




observables, mainly in [l| and decay rate j^, in —)■ decay rate [3| and in the 

ratio Rk = BR(R —)■ iP/i+/r“)/BR(R —)■ Ke^e~) jd]. Interestingly all these deviations are 
associated with the quark level transition b —?• . 

In this paper, we would like to focus on the semileptonic decay mode R° —)■ )■ 

R7r)/i+/x“ which is quite an interesting channel, as the measurement of four-body angu¬ 
lar distribution provides a large number of observables which can be used to probe and 
discriminate different scenarios of NP. Theoretical predictions for such observables are par¬ 
ticularly precise and free from hadronic uncertainties in the low-range of dimuon invariant 
mass squared i.e., 1 < < 6 GeV^. While the observed forward-backward asymme¬ 

try is systematically below the SM prediction, the zero crossing point is consistent with it. 
Also there are few other deviations from the SM expectations have been observed by LHCb 
experiment in the angular observables. The largest discrepancy of 3.7cr encountered in the 


observable [l| in the bin g^ G [4.3, 8.68]. Another interesting observable to look for new 
physics is the isospin a^mmetry distribution, which is measured by LHCb experiment in 
the entire g^ spectrum |5| . The leading uncertainties in the R —)■ K* form factor is expected 
to cancel in this asymmetry. 

The angular distributions of R —)■ K*Pl~ processes with the dilepton invariant mass 
has been studied by various experiments such as BaBar, Belle, CDF, and LHCb. All these 
experiments cover the full kinematical dilepton mass region, i.e. Amf ^ g^ ^ 
leaving the regions around g^ ~ general the kinematically allowed region 
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can be classified into three regions and different theoretical approaches nsnally adopted to 
stndy the properties of different observables. In the region of large hadron recoil i.e., for 
^ the kaon is very energetic and varions physical observables can be compnted 

nsing QCD factorization (QCDf) approach. The intermediate valnes of i.e. 7 GeV^ ^ 
^ 14 GeV^ fall into the narrow-resonance region and cuts are employed to remove the 
dominated charmonium resonance (cc) = J/V’, 4’’ backgrounds from B —)• K* (cc) K*l^l~. 
The larger dilepton invariant mass region i.e., q^ ^ 14 GeV^ corresponds to low-recoil limit 
and in this region the kaon energy is around a GeV or below. Here soft collinear effective 
theory (SGET) and QGD factorization approaches are not justified properly and become 
invalid near the zero recoil point = {jnB — rriK*)'^- At low recoil the heavy to light 


decays can be studied by an operator product expansions in 1/Q where Q = ymi,, yq^j 
i.e. is of the order of the mass of the b quark, rrib a, l7|. The combination of operator 
product expansion (OPE) with the heavy quark effective theory (HQET) and the use of 
improved Isgur-Wise form factor relations 6|, ISj allows to obtain the B -y K*l^l~ matrix 
element expansion in the strong coupling and in power corrections suppressed by the heavy 
quark mass, in low recoil. 

nr 

Recently, the observed anomalies associated with b —)■ sl~^l processes at LHGb [lH4| 
have attracted a lot of attention to look for new physics both in the context of various new 


physics models as well as in model-independent ways |9|-[l3|. In this paper, using scalar 
leptoquark model, we would like to study the B^ —)■ processes, which contain quite 

a large number of clean observables in the full kinematics except the intermediate q^ region. 
In particular, we are interested to look for the effect of scalar leptoquark on some of the 
observables such as dilepton mass spectra, lepton-angle distribution and various asymmetries 
like forward-backward asymmetry and isospin asymmetry. 

The similarities between leptons and quarks lead to the fact that there could exist lep- 
toquarks (LQs), which are color triplet bosons and carry both lepton (L) and baryon {B) 
quantum numbers. Leptoquarks violating both B and L numbers are generally considered 
to be very heavy at the level of 0(10^®) GeV to avoid proton decay. On the other hand LQs 
conserving B and L can be light and can have implications in the low energy phenomena. 


The existence of leptoquar 


is has been proposed in many extensions of the SM e.g.. Grand 


Unified T 
scenarios 


leories (GUTs) 


14j |. Pati-Salam model |^, tec hn icolor models 


Q. 


composite 


171], etc. The spin of leptoquarks could be either one (vector leptoquarks) or zero 
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(scalar leptoquarks). Scalar leptoquarks are encountered in extended technicolor models 
and models with compositeness of quark and lepton [^, 171 at TeV scale. However, in this 
case the bounds from proton decays may not be relevant and leptoquarks may give signa¬ 


tures in other low energy processes 


18| . The phenomenology of scalar leptoquark and the 


contribution to new physics has been quite well studied in the literature 


18 


22| . However, 


the effect of scalar letoquarks in various observables associated with B —)■ process 


is not yet explicitly studied. In Ref. 


20| model independent constraints on leptoquarks from 


b —>■ processes are obtained. In this paper, we would like to see how the scalar lep¬ 

toquarks affect these observables and whether it would be possible to differentiate between 
these two scalar LQ models from some of these observables. 

The plan of the paper is as follows. We present a brief discussion on the effective 
Hamiltonian for b —)■ processes in the SM as well as in leptoquark model in Section H. 

The new physics contributions to these processes due to the exchange of scalar leptoquarks 
and the constraint on leptoquark parameter space from the rare decay mode Bg —)■ 
have also been discussed. The constraints obtained from Bg — Bg mixing is discussed in 
Section HI. The observables associated with the decay modes B^ —)■ are presented 

in Section IV. Our predicted results on branching ratio, isospin asymmetry parameter and 
various form factor independent observables in the angular distribution are also presented 
in this section. Section V contains the summary and conclusion. 


II. EFFECTIVE HAMILTONIAN FOR b sl+B PROCESS 


The effective Hamiltonian describing the flavour-changing quark level transitions b —)• 
in the standard model is given as 231 

VtbVtg 

i=l 


^e// 


x/2 


b 

(s(T^y{mgPL + mbPR)b^F^'' 


+Cl"^{srPLb)h,l + C,o^{srPLb)h,75l 


a 


( 1 ) 


where Vqqi denote the CKM matrix elements, Gp is the Fermi constant, a is the hne- 
structure constant, Pl,r = (1 T 75)/2 is the chirality projection operator and Q’s are the 
Wilson coefficients. The values of the Wilson coefficients evaluated at the scale /i = rrib in 
the next-to-next-leading order are listed in Table-1. 
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TABLE I: The SM Wilson coefficients evaluated at the scale /i = 4.6 GeV 



G2 

Cs 

G4 

Cs 

Ce 


r^ff 

C^9 

Cio 

-0.3001 

1.008 

-0.0047 

-0.0827 

0.0003 

0.0009 

-0.2969 

-0.1642 

4.2607 

-4.2453 


The effective Hamiltonian described above in Eq. ([T]) will receive additional contributions 
arising due to the exchange of leptoquarks. We will present the modihed Hamiltonian in 
the presence of scalar leptoquarks in the subsection below. 


A. New physics contribution from scalar leptoquark 


Models with scalar leptoquarks can modify the effective Hamiltonian due to the exchange 
of leptoquarks and will give measurable deviations from the predictions of the SM in the 
favor sector. Here we will consider the minimal renormalizable scalar leptoquark model 

n 

18l |. containing one single additional representation of SU{3) x SU{2) x U{1) which does 
not allow baryon number violation in perturbation theory. There are only two such models 
which are represented as A = (3, 2, 7/6) and X = (3, 2,1/6) under the SU{3) x SU{2) x U{1) 
gauge group. Here, we are interested to study the effects of these scalar leptoquarks which 
potentially contribute to the quark level transition b —)■ sl^l~ and constrain the underlying 
couplings from experimental data on Bg —)■ Although the details of this method has 
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26|, here we will briefly mention about the main points for 


been discussed in Refs, 
completeness. 

The interaction Lagrangian for the scalar leptoquark X = (3, 2, 7/6) couplings to the 


fermion bilinear 


18| is 


C = - A«4.YtQi + h.c. , 


( 2 ) 


where i,j are the generation indices, Ql {Li) is the left handed quark (lepton) doublet, X 
is the scalar leptoquark doublet, ur (cr) is the right handed up-type quark (charged lepton) 
singlet and e = ia 2 is a 2 x 2 matrix. 

After expanding the SU{2) indices and performing Fierz transformation, the contribution 
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to the interaction Hamiltonian for the process b —)■ is 


^32^22* ^32^22* 

'^lq = g^2 - 75)&] [h7/.(l + 75)h] = (^9 + Oio), (3) 


Ly -tiv±Y 

which can be written analogous to the SM effective Hamiltonian as 

Hio = (cro,+ . 

y I'K 

Thus, one obtains the new Wilson coefficients 


(4) 


riNP _ r^NP _ 
*-^9 ~ '^10 ~ 


TT 


\32 \ 22 * 


(5) 


2V2GfaVtbVZ M2 

Similarly, the corresponding Lagrangian for the coupling of scalar leptoquark X = (3, 2,1/6) 


to the fermion bilinear is 


C = + h.c.. 


( 6 ) 


Proceeding in the similar manner as done in the previous case, the interaction Lagrangian 
becomes 

\ 22 \ 32 * \ 22 \ 32 * 

y-LQ = [57^(1 + 75)&] [h7.(l - 75)h] = (C?9 - C?'io), (7) 

where O'^ and O'^q are dimension-six operators obtained from Og and Oiq by the replacement 
Pl -f-)- Pr and their respective new Wilson coefficients due to the exchange of the leptoquark 
X = (3, 2,1/6) are given as 


y-ilNP _ yy/NP _ 
'-'9 ~ '-'10 ~ 


TT 


^22^32* 


( 8 ) 


2V2GfaVttV,l M2 

After having the new Wilson coefficients in hand, we now proceed to constrain the combi¬ 
nation of LQ couplings by comparing the theoretical 27|] and experimental branching ratios 


28 


30| of n , as these new coefficients contribute to the Bg —)■ ijG/ i process as 


well. Furthermore, we require that each individual leptoquark contribution to the branch¬ 
ing ratio does not exceed the experimental result. The constraint on leptoquark parameter 


space has been extracted in 
The allowed region in r 
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26| . therefore here we will simply quote the results. 

plane which is compatible with the Icr range of the 
experimental data is 0 < r < 0.1 for the entire range of , i.e.. 


0 < r < 0.1 , for 0 < 


lNP 


< 2tt 


(9) 
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where r and are defined as 


- C[^P)/C- 
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( 10 ) 


However, in onr analysis we will nse relatively mild constraint, consistent with both mea- 
snrement of BR(i?s —)■ and BR(S2 —)■ XsfiPfi~) [2^ as 


0 < r < 0.35 , with 7r/2 <0 < 3tt/2 . 


( 11 ) 


It shonld be noted that the nse of this limited range of CP phase, i.e., {n/2 < < 37r/2) 

is an assnmption to have a relatively larger valne of r. The constraint on r can be translated 
to obtain the bonnds for the leptoqnark conplings using Eqs. (5), (8) and (11) as 

1 ^ 32 ^ 22*1 


0 < 




< 5 X 10 ^ GeV for vr/2<0 <37r/2. 


( 12 ) 


III. BOUND FROM B, - B, MIXING 


Now we will obtain the constraint on the leptoqnark couplings from the mass difference 
between the meson mass eigenstates (AMg), which characterizes the Bg — Bg mixing 
phenomena. In the SM, Bg — Bg mixing proceeds to an excellent approximation through the 
box diagram with internal top quark and W boson exchange, and the effective Hamiltonian 
describing the AB = 2 transition is given by 311 




lOvr^ 


X; M^So{xt)r]B{sb)v-A{sb)v-A , 


where Xt = rjB is the QCD correction factor and So{xt) is the loop function 

4:Xt — llacj + Xj 3 logxtx^ 


So{xt) = 


4(1 - xty 


2il-x,f ’ 


(13) 


(14) 


with Xt = Thus, the Bg — Bg mixing amplitude in the SM can be written as 

where the vacuum insertion method has been used to evaluate the matrix element 
(aKsi")! - 75)6)(S7„(1 - 75)4)|S.) = ■ 


(15) 


(16) 


The corresponding mass difference is related to the mixing amplitude through AMg = 
2|Mi2|. Now using the particle masses from [32|, tjb = 0.551, the Bag parameter ISb^ = 
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1.320 ±_0.017 ± 0.030 the decay constant fss = 225.6 ± 1.1 ± 5.4, t-quark mass rrit = 165.95 
from 


33| . we obtain the value of AMg in the SM as 

= (17.426 ± 1.057) ps-\ 


which is in good agreement with the experimental result 321 

AM, = 17.761 ±0.022 ps"^ 


(17) 


(18) 


However, the central value of the theoretical prediction deviates from the corresponding 
experimental value. The ratio of these two results yields 

AM,/AMf^ = 1.019 ± 0.062 , (19) 


which is consistent with one, but it does not completely rule out the possibility of new 
physics in Bg — Bg mixing. The mixing amplitude receives additional contribution due to 
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FIG. 1: Box diagram for Bg — Bg mixing phenomenon in the leptoquark model. 


the flow of leptoquark and charged lepton/neutrino in the box diagram as shown in Fig.l. 
For X(3, 2, 7/6) LQ, there will be contribution coming only from charged lepton in the loop 
whereas for X(3,2,1/6) both charged lepton and neutrino will contribute to the mixing 
amplitude. 

The effective Hamiltonian due to the leptoquark X(3,2,7/6) and charged lepton in the 
loop is given by 


(yk^si*)2 ^ I' ^2 


= E > (i|j . 

where the loop function I{x) is given as 

1 — + 2x log X 


( 20 ) 


I(x) = 


(1 — xy 


( 21 ) 
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which is always very close to J(0) = 1. For X(3,2,1/6) contribution there will be charged 
lepton as well as neutrinos in the loop and the corresponding effective Hamiltonian becomes 




1287r2 


2=e,/i,r 


Ml 


rrij 

Ml 


+ 


Mh 


{brPRs){h^.PRs) . 


( 22 ) 



FIG. 2: The allowed parameter space for the leptoquark couplings in the jA^^A^^/M^I vs. (j)^^ 
plane obtained from the mass difference between Bg meson mass eigenstates {AMg). The left panel 
corresponds to bounds on X(3, 2,7/6) and right panel is for X(3, 2,1/6) couplings. 


To obtain the constraints on the leptoquark coupling, we require that individual lepto¬ 
quark contribution to the mass difference does not exceed the la range of the experimental 
value. Since we are interested to obtain the bounds on A^^ and A'^^ couplings, we consider 
the muon contribution to the mixing amplitude. Neglecting the mass of muon and using 
Eq. ffT6l) . we obtain the contribution due to leptoquark exchange as 

/ \ 32* \ 22A2 

Ml? = 'IbBbJIMb, , for X(3,2,1/6) 

Ml? ^ , forX(3,2,7/6). (23) 

Thus, including both SM and leptoquark couplings the total contribution to mass difference 
is given as 


AMg = AM: 


SM 


C /(A32A22*)2 

^ 16GlV,lV,fml,S„(x,} V Ml 


(24) 


where the constant c = 1 for X(3, 2,1/6) and 1/2 for X(3, 2, 7/6). Now varying the ratio of 
mass difference {AMg/AM^^) within its la allowed range flT^ . we obtain the constraint on 


9 


















/Ms\ as shown in Fig. 2, where the left plot corresponds to constraint on X(3, 2, 7/6) 
and the right plot shows the constraint on X(3,2,1/6) conplings. From the hgnre, the 
bonnds on \X^‘^X^‘^/Ms\ for the entire allowed range of (j)^^ are fonnd to be 


0 < 

0 < 


^ 32^22 


Ms 

^ 32^22 


< 7.5 X 10-® GeV 


T-l 


Me 


< 5.0 X 10-® GeV 


T-l 


for X(3,2,7/6), 
for X(3,2,l/6). 


(25) 


It should be noted that using the Bg — Bg mass difference, we obtained the bounds on 
lA^^A^^I/Mg, whereas using the Bg —)■ and B^ —>■ X^/x/i data the bounds on |A^^A^^|/M| 


have been obtained. So to correlate these two resu 
scalar leptoquark Ms- Recently GMS collaboration 


ts, we need to know the mass of the 


34| with 8 TeV data set excluded the 


hrst generation leptoquarks with masses less than 1010 (850) GeV for /3 = 1.0 (0.5), where 
(3 is the branching fraction of a leptoquark decaying to a charged lepton and a quark. 
The second generation scalar LQs are excluded with masses less than 1080 (760) GeV for 
(3 = 1.0 (0.5). They also ruled out at 95% conhdence level the single production of hrst 


generation LQs with coupling and branching fraction o 
and for second generation for masses below 530 GeV 35 


1.0 for masses below 1730 GeV 
. ATLAS Gollaboration excluded 


at 95% G.L. the scalar leptoquarks with masses upto 1050 GeV for hrst generation LQs 
(i.e., rriLQi < 1050 GeV), mr, Q 9 . < 1000 GeV for second generation LQs and < 625 
GeV for third generation LQs 36|. Hence, if we scale the couplings obtained from Bg — Bg 
mass diherence for a benchmark leptoquark mass of 1 TeV, the bounds in Eq. fl25ll can be 
translated as 


0 < 

0 < 


A32A22 


Ml 

A32A22 


< 7.5 X 10"® GeV-^ 


Ml 


< 5.0 X 10"^ GeV 


-2 


for X(3,2,7/6) 
for X(3,2,l/6). 


(26) 


Since these bounds are reasonably higher than those of obtained from Bg —)■ fifi and B^ —)■ 
Xgfifi , we will use the bounds flT^ as mentioned in the previous section, in our analysis. 


IV. ANALYSIS OF K*l+l- PROCESSES 

Here we will consider the decay modes —>■ K*l^l~. At the quark level, these processes 

proceed through the FGNG transition b —)■ which occurs only through loops in the 
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SM, and hence they are quite suitable to look for new physics. Moreover, the dileptons 
present in these processes allow one to formulate several useful observables which can serve 
as a testing ground to decipher the presence of new physics. 

The transition amplitude for these processes can be obtained using the effective Hamil¬ 
tonian presented in Eq. ([T]). The matrix elements of the various hadronic currents between 
the initial B meson and the hnal K* vector meson can be parameterized in terms of seven 


form factors by means of a narrow width approximation. The relevant form factors 
given as 


371 are 


{K* (pk*) \s'yf,PL,Rb\B (p)) = 

-ie*-q){ 2 p- q)^ 


V{s) 


rriB + rriK^ 
y42(s) 2mK^‘ 


=F X e* {mB + mK*)Ai{s) 


rriB + TriK* s 
{K* (pk*) \Ma^,yq''PL,Rh\B (p)) = -ie^yafse‘'*p°‘q^ Ti{s) ± 


e*-q) [^ 3 ( 5 ) -Ao{s)]q^ |, 
1 




-{e*-q){2p-q)f, T 2 {s) + {e*-q) 


% 


(m| -m^*) 


( 2 p - q)^ 


rriK*) 


^ 3 ( 5 ), (27) 


where q = pi+ +Pi-, s = q^ and is the polarization vector of K*. The form factors Hi, 
H 2 and H 3 are related to each other through 


{rriB + mK*) (m^ 

^ 3 ( 5 ) =-- 


UIr* 




(28) 


2 mK* 2 mK* 

The amplitude for the process B —)■ K*{-^ K'k)PI~ can be represented by seven transversity 
amplitudes, and The explicit form of these amplitudes (up to corrections 0{as) 

are presented in Appendix A (B) for low q^ (high region. 

Assuming the K*^ K~'k^ to be on the mass shell, the full angular distribution of the 

decay B —)■ K*^ (—)■ Pl~ can be described by four independent kinematic variables, 

the lepton-pair invariant mass and the three angles 6 '^*,_^and 0. The differential decay 


distribution in terms of three variables can be written as 

d^r 9 




(29) 


dq^ d cos 61 d cos 6 k* d(j) 2>2 ti 
where the lepton spins have been summed over. Here q^ is the dilepton invariant mass 
squared, 61 is dehned as the angle between the negatively charged lepton and the B in the 
dilepton frame. Or* is the angle between K~ and B in the K~ti^ center of mass system and 
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0 is given by the angle between the normals of the K and the dilepton ) planes. 
The full kinematically physical region phase space is given by 

4mf ^ ^ {ruB — rriK*)^ , —1 ^ cos6*/ ^1, —1 ^ cosOk* ^1, 0^0^ 27r. (30) 

where rriK*, fni are the masses of B meson, K* and lepton respectively. More explicitly 

the dependence of the decay distribution on the three angles can be written as 

J (g^, 6 u 6 k* , 0 ) = <^1 sin^ 6k* + Ji cos^ 6k* + {J 2 sin^ 6k* + J 2 cos^ 6k*) cos 26 i 

+ J 3 sin^ 6k* sin^ 61 cos 20 + J 4 sin 26 k* sin 26i cos 0 + J 5 sin 26 k* sin 61 cos 0 

+ ( Jg sin^ 6k* + Jq cos^ 6k*) cos 61 + J 7 sin 26 k* sin 61 sin 0 
+ Jg sin 26*x* sin 26*/sin 0 + Jg sin^ sin^ 0/sin 20 , (31) 

where the coefficients = 0*"“^ (g^) for i = 1,....,9 and a = s,c are functions of the 
dilepton mass, and are expressed in terms of the transversity amplitudes Aq, Aj_, and 
At as given in Appendix C. 

The dilepton invariant mass spectrum for B —)■ decay after integration over all 

angles (38| is given by 

where Ji = 2J^ + J^. An interesting observable to look for new physics is the zero crossing of 
forward-backward asymmetry Afb, which can be obtained after integrating the 4-differential 


distribution over 0 and 6k* angles and is dehned as 
-^FB {<f) = 


r , n f 

IA 

3 Jq 


acos&i- 


' dq'^d cos 61 


>dT 

dq^ 


(33) 


8 dV/dq^ ■ 

The longitudinal and transverse polarization fraction of the K* meson can be dehned in 
terms of the transversity amplitudes as 


Fris) = 




I^Op + I^IlP + \^± 


2 ’ 


Fris) = 


l^llP + I^J 


l^oP + l^iiP + \^± 


2 ’ 


(34) 


and in terms of the angular coefficients JJs these observables can be expressed as [41 1 

3 Jf - J| 


Fl (s) = 


AdV/dq^ ’ 


Ft{s) = 1 - Fl{s) , 


(35) 
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so that one can define the ratio of K* polarization fraction ok* as ^ 

-J2 


ax* (s) = -1 = 


VI 


( 36 ) 


The transverse asymmetries are given as 42 1 

yiy(s) = 


4?' (s) = 


a!^ (s) = 


IVP + 


J3 


2J|’ 


A(jV 

+ f^i iJr)'^ 

V - 2 J| 

(2J| + J3) 

(HJs)^ 

+ f^i iJ^f 

V 

to 

+ Pf iJr)'^ 

\A^Af* 


IVP 

+ l^llP ’ 

J9 



1/2 


1/2 


aP (s) = 

i^±r i^iir 

Another set of interesting observables are the six form factor independent (FFI) observables 
431], which are given by 

J 3 


Pi is) = 
P,is) = 
Pe is) = 


2 J| ’ 


V2Ji 


P2M 

7 ® 

- R ^ 



(s) = 

I 3 iJ 3 

J3) ' 

2 J| ( 2 J| + J3) 




y-2j|(2j|-J3)' 


(38) 


A slightly modified set of clean observables -P 4 5 g which are related to Pa, 5,6 through the 
relations {q] 


Pi = PaV^ - Pi = 
F' = FsV + ^i = 


Ja 


J 3 


F' = F6 V-^i = 


2 \/ J 2 J 2 

—Jr 


(39) 


A. Observables in the large recoil 

After getting familiar with the different observables, we now proceed to study these 
observables in the large recoil limit. For that we need to know the associated form factors 
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for B —)• K*l^l~ process. In the heavy quark limit the QCDf form factors obey symmetry 
relations and at leading order in the 1/E expansion, they can be expressed in terms of two 
universal soft non-perturbative form factors and In order to calculate the universal 


form factors we use the QCDf scheme 
tub 
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niB + rriK* 


V{q^ 


44l |. where they are expressed as 

niB — rriK* 


C (TP \ + "mK* , , 2\ 

i\\{EK*) = - Ai{q ) 


2E 


K* 


rriB 


A2{q^). (40) 


The dependence of the form factors V, Ai, A 2 can be parameterized as 


= 

= 

^2 (g^) = 


ri 


q‘^/m\ 

r2 


+ 


r2 




q^/m).^ 

ri 


+ 


r2 




(41) 


The values of the parameters involved in the calculation of form factors are taken from 45 |. 
The Tj formfactors are related to the universal form factors as 


2E 


T,{q^) = UEk*) , T2{q^) = ^^±{Ek*) , Uq^) = UEk*) - ^ii{Ek*) • (42) 


tub 


TABLE II: Summary of the values of various input parameters. 


a 

1/137 

0-1,± 

0.10 

Ois{Mz) 

0.1184 

0-2, ± 

0.13 


4.6 GeV 

«l,ll 

0.10 

'^c,pole 

1.4 GeV 

02,11 

0.09 

fK*,± 

0.185 GeV 

A 

0.22537±0.0006 

fK*,\\ 

0.220 GeV 

A 

n SI /1+0-023 

U.5i4_o 024 

fs 

0.2 GeV 

P 

0.117 ± 0.021 



P 

0.353 ± 0.013 


After getting familiar with the different observables and the associated form factors for 
B —)■ K*l^l~ processes in the high recoil limit, we now proceed for numerical estimation. 


The masses of particles and the lifetime of B meson are taken from 


32l | . For the leptoquark 
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couplings we use a representative value for the parameter r as r = 0.3 and vary the associated 
phase between 7r/2 < cf)^^ < 37r/2. Furthermore, we will present most of the the results 
only for X(3, 2, 7/6) LQ and only a few representative plots for X(3, 2,1/6). The values of 
quark masses and all the input parameters used in our analysis are listed in Table-II. 




FIG. 3: The variation of branching ratios of i? —?■ (left panel) and B —?■ K*e'^e (right 


panel) with low in the standard model and in leptoquark mode. The green bands are due to the 


leptoquark contributions and the gray bands represent the theoretical uncertainties from the input 


parameters in the SM. The plots in the top panel are for X(3,2,7/6) and those in the bottom 
panel are for X(3,2,1/6). 


In Fig. 3, we show the variation of the branching ratios for B —)■ (left panel) 

and B —)• K*e^e~ (right panel) in the low region. The plots in the top panel are 
for X(3,2,7/6) and those in the bottom panel are for X(3,2,l/6). The variation of the 
longitudinal and transverse polarization fractions of K* has been shown in Fig. 4 and that 
of forward-backward asymmetry in Fig. 5. From these hgures one can see that the affect 
of the LQs X(3, 2, 7/6) and X(3, 2,1/6) are quite different and one can easily differentiate 
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between these two models from the measured values of the K* polarization fractions FL{q^) 
and The transverse asymmetry parameters A^\ A'^\ Aim and K* polarization factor 

ax* variations with are presented in Fig. 6. The variation of form factor independent 
observables as a function of dimuon invariant mass squared have shown in Fig. 7. The 
total branching ratios and the asymmetries integrated over the range G [1,6] GeV^ are 
presented in Table III and the allowed range of transverse asymmetry and the form factor 
independent observables are given in Table IV. It should be noted that in the leptoquark 
model these observables deviate signihcantly from their SM predictions. 




FIG. 4: Variation of longitudinal and transverse K* polarization fractions in the low (f region. 
The plots in the top panel are for X{3, 2, 7/6) and those in the bottom panel are for X{3, 2,1/6). 


Another interesting observable is the Isospin asymmetry distribution, which has been 
recently measured by LHCb experiment at 3 fb“^ data set This asymmetry arises due to 
the non-factorizable part where photon is radiated from the spectator quark in annihilation 
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FIG. 5: The variation of forward backward asymmetry for X(3,2,7/6) (left panel) and for 
X(3,2,l/6) (right panel) with q^. 

TABLE III: The predicted values for the integrated branching ratio, forward-backward asymmetry, 
isospin asymmetry and the polarisation fractions in the range G [1,6] GeV^ for the B —)• 
process. 


Observables 

SM prediction 

Values in LQ model 

A(3,2,7/6) 

Values in LQ model 

A(3,2,l/6) 

BR(B ^ 

) (7.738 ± 0.464) x 10“ 

^ (6.88 ^ 8.73) X 10-7 

(4.8 ^ 8.3) X 10-7 

BR{B K*e+e-) 

(7.742 ± 0.465) x 10“ 

7 (7.03 ^ 8.73) X 10-7 

(4.85 ^ 8.38) X 10-7 

(Afb) 

-(0.09 ±0.005) 

-0.11 ^ 0.004 

-0.185 ^ -0.08 

(Ai) 

-(0.03 ±0.002) 

-0.06 ^ -0.04 

-0.02 ^ -0.01 

(Fl) 

0.71 ± 0.043 

0.7^ 0.8 

0.5 ^ 0.7 

{Ft) 

0.29 ± 0.017 

0.2 ^ 0.3 

0.3 ^ 0.5 


and spectator scattering. The CP-averaged isospin asymmetry is defined as 



dAi _ dr [B^ K*H+l-] /dq^ - dr [B^ iC*±/+/-] /dq^ 

~df ~ dr [B^ K*H+l-] /dq^ + dr [5± ^ K*^l+l-] /dq^ ' 


(43) 


Including longitudinal photon polarizations appearing for q"^ ^ 0, the isospin asymmetry 
distribution in the QCD factorization scheme is given by 


Aj{q^)=Re{bj{q^)-bi{q^)) 


(g2) 


_ 

|Gio()U,)|2 + |Cf^(g2)|2 G(g2 


(44) 
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with 


F{q^)=l + 


and 


G (g2) = 1 + 


1 El,ml eg jq^) {Q)-bi{q)) [C'^ll (^2) |2 

4 q^m\. ^2 (^2) ^^2)) |^W± ^^2) |2 

1 (?^) 1^9°^" (g^) 1^ + \Cio (/if,) p 


4 q^ml. (g2) (^2) |2 + |2 

where the generalized standard model Wilson coefficients are 

2mbmB 


CE (9=) = c, (W) + r if) + M 


9 


and 


d“’" (i-) = C, (Mt) + F («^) + iii,) 


The 6“ (a =T, ||) terms appearing in the above eqnations are given as 


i-t if) = 


2ATT‘^mBfBeq 

q^u (9=) c'E (9") V ™» 




/q (,= 


6 Ab,+ (g^) rriB 1 — g^/m' 


B 


and 


i>! (9=) 


2 Ati‘^ fBCgTriK- 


fx* 


E if) 


(45) 


(46) 


(47) 


(48) 


(49) 


(50) 


( 9 =) Cf" (g 2 ) V3Ab,- ( 9 ") 

where the expressions for the terms K ^2 presented in Appendix E. The variation of 
isospin asymmetry distribntion with respect to dimnon invariant mass sqnared has given in 
right panel of Fig. 8 and Table III contains the allowed range of isospin asymmetries. 


B. Observables in the low recoil 

At low recoil the exclnsive B decays depend on the improved form factor 

relations and an operator prodnct expansion (OPE) in 1/Q. The OPE controls the non- 
pertnrbative contribntions from fonr-qnark operators and is important for charm qnark, 


whose operators enter withont any snppression from CKM matrix elements and sma! 
coefficients. The QCD operator identity for massless strange qnark (m^ = 0) is [6 

id'' {sia^yh) = -nib (s'liib) + id^ (sb) - 2 , 


1 Wilson 


49| 


(51) 


which allows to extract relation between the form factors Ti and V and the matrix elements 
of the cnrrent si^^b. The improved Isgnr-Wise relations to leading order in l/nib inclnding 
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TABLE IV; The predicted integrated values of the FFI observables and the CP violating observables 
in the range G [1,6] GeV^ for the B —)■ process. 


Observables 

SM prediction 

Values in LQ model 

V(3,2,7/6) 

Values in LQ model 

V(3,2,l/6) 

{Pi) 

-0.044 ±0.003 

-0.037 ^ -0.046 

-0.017 ^ 0.14 

{P2) 

0.203 ± 0.012 

0.08 ^ 0.23 

0.19 ^ 0.21 

(Ps) 

-(6.0 ±0.4) X 10“^ 

-(9.4 ^ 1.8) X 10-3 

-0.014 ^ 0.08 

{Pa) 

0.395 ± 0.024 

0.294 ^ 0.45 

0.24 ^ 0.52 

{P5) 

-0.204 ±0.012 

-0.42 ^ -0.13 

-0.39 ^ -0.13 

{Pe) 

-0.0075 ± 0.0005 

-0.07 ^ 0.08 

0.079 ^ 0.112 

(4b 

0.55 ± 0.033 

0.56 ^ 0.6 

0.37 ^ 0.6 

<4'> 

0.87 ± 0.05 

0.82 ^ 0.94 

0.99 ^ 1.5 

{Aim) 

(1.7 ±0.1) X 10“^ 

(2.0 ^ 5.3) X 10-^ 

-0.023 ^ 0.005 

{otK*) 

3.73 ± 0.22 

3.6 ^ 5.6 

1.8 ^ 3.6 


radiative corrections are 


Ti (q^) = nV (g2) , 

T2 (q^) = nA, (g^) , 

Ts (g^) = .A, (g^) ^ , (52) 


where 


(/^) (/i) 


K = 1 + 


G 


(.A 


0 (h) 


niB 


(53) 


and at /i = mb inclnding 0(as) corrections, it reads k = 1 + 0{al). We have shown the 
variation of branching ratio for B —)■ (left panel) and B —)■ K*e^e~ (right panel) 

with g^ in Fig. 9. In the low recoil region the variation of forward-backward asymmetry, 
isospin asymmetry, longitndinal and transverse polarization fractions of K* with respect to 
g^ are shown in Fig. 10. Fig. 11 shows the variation of P 2 , with respect to dimnon 
invariant mass squared. Table V contains the integrated values of branching ratio, forward- 
backward asymmetry and isospin asymmetry in the low recoil i.e. q^ G [14.2,19] GeV^. It 
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FIG. 6; The variation of transverse asymmetry parameters , Ai^ and the K* polarization 

factor with in high recoil for X(3,2,7/6) LQ. 


should be noted from these hgures that at high there is no signihcant deviation between 
the SM results and the leptoquark predictions. 


V. CONCLUSION 

In this paper we have studied the rare semileptonic decays —)• using the 

simple re-normalizable leptoquark model in which a single scalar leptoquark is added 
to the standard model with the requirement that proton decay would not be induced 
in perturbation theory. The leptoquark parameter space is constrained using the recent 
measurement on Bg —)■ Using such parameter space we obtained the bounds on 

the product of leptoquark couplings. We then estimated the branching ratios, isospin 
asymmetries and forward-backward asymmetries for B —)■ process in full physical 

region except the intermediate region of q^. The CP violating observables and the form 
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q^[GeV^] 



FIG. 7: The variation of the observables 77,2,3 and with for X(3, 2,7/6) LQ. 


factor independent observables have also been studied in the leptoquark model. It is found 
that in these models, there could be significant deviations in these observables in high recoil 
but comparatively less in low recoil regimes. We found that the time-integrated values 
of some of the asymmetry parameter have deviated signihcantly from their corresponding 
SM values, the observation of which in the LHCb experiment would provide the possible 
existence of leptoquarks. 
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TABLE V: The predicted values for the integrated branching ratio, forward-backward asymmetry 
and isospin asymmetry in the range E [14.2,19] GeV^ for the B —)• K* process. 


Observables 

SM prediction 

Values in LQ model X{3, 2, 7/6) 

BR{B AV+p-) 

(8.5 ± 0.51) X 10-^ 

(8.93 ^ 9.31) X 10-^ 

BR{B K*e+e-) 

(8.52 ±0.511) X 10"'^ 

(8.92 ^ 9.32) X lO-"^ 

{Afb) 

0.4 

0.34 ^ 0.38 

{Ai) 

-(2.75 ±0.17) X 10-3 

(-3.3 ^ 2.5) X 10-3 

{P 2 ) 

-0.42 

-0.41 ^ -0.38 

(4*') 

0.57 

0.53 ^ 0.58 

(Fl) 

0.35 

0.35 

(Ft) 

0.65 

0.65 ^ 0.66 




123456 123456 

q^[GeV^] q^[GeV2] 


FIG. 8: The variation of isospin asymmetry for A(3, 2, 7/6) (left panel) and for X(3, 2,1/6) (right 
panel) with q^. 
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FIG. 9: The variation of branching ratio for B —>• with high for X(3, 2, 7/6) LQ. 






FIG. 10: The variation of forward backward asymmetry, isospin asymmetry, longitudinal and 
transverse polarization fraction with high (f‘ for X(3,2,7/6) LQ. 
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FIG. 11: Variation of the observables P 2 ) and for V(3,2,7/6) LQ. 

Appendix A: Transversity amplitudes at NLO in the large recoil 


In the large recoil limit the transversity amplitudes at next to leading order (NLO) within 
QCDf can be given as 38, 

= iV\/^ 


((C, + + C'«-) T (C„ + c«r + CID + ?!^r, 


+ mK* q 


2 '-L 


= — A\/2 — m^*) 


((Gs + c,'"’ - CfT (Cio + C,y - c:r)) 


Ai 


niB — mx* 


Arub Ex* 


T; 


aL,R _ 
/Iq — 


X 


N 


tub 

((09 + - of T (Olo + o(^^ - o(r)) 

(f) 


2mx*\f^ 

(m| - - q^) {jub + mx*) Ai [q^] - A 


rriB + mK* 


2mb 


2E 


K* 


rriB 


{rn\ + - q^) 7] 


A 


J_ 2 2 


rr + 'O- 


/i, = iVihiOv/A (c,„ + c,y - c;d , 

f g2 mic* A[| 


(Al) 


A5 = 0, 

where Of^ and Of q are the new Wilson coefficients arising due to leptoquark exchange 
and E^ is the energy of the kaon in the B meson rest frame and is given as 

(A2) 


Ex* — 


m% + — q'^ 


The normalization constant N is given as 


N = 


Gl 


a 


3 . 2 ^^'K^mB 


2m B 


■\VtbV:fs^^\|3l 


1/2 


(A3) 
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where 


A = m% + — 2 + rn\q^^ , and s = g^/ 


m 


B ■ 


(A4) 


The transversity amplitude At contains Ay and negligible for massless lepton. It contributes 
only for mi 7 ^ 0. The light-cone distribution amplitude for K* is given by 

= 6 n (1 - u) {1 + ai {2u - 1) + as {R*){2u - 1 )} , (A5) 


where the moments are 


Jo ^ 




o-Al{mBUJo) 


Uq 


[-Ei (g^/msWo) + ivr] . 


(A 6 ) 


(A7) 


The detailed expression for the function 7^ (a =T, ||) at NLO in the QCDf framework is 
given in Appendix D. 


Appendix B: Transversity amplitudes in the low recoil 

The transversity amplitudes to leading order in 1/mb at low recoil are given as 

2 m, 


A L,R 

A_^ = i 


{(Cfj> + cr + cr) =F (C.„ + Cf/ + CiD) + 


/±. 


\ LiR 

An' = —I 


A 

Aq’ = 


(C'l" + T (C'lo + - C'JD) + 

(C‘" + cr - Cfn =F (C.o + - c;r)) + 


f\\ ’ 

/o, (Bl) 


where the form factors read 


/o = NruB 
and the normalization factor is 


= Nms V, 

1 -|- mx* 

/|l = NmB\f2 (1 -1- niK-) Ai, 

{1 - s - ml-*) (1 + ttikA Ai - AAs 


N = 


2mK* (1 + Aik*) Vs 

n^~\ 1/2 

Gpa'^\Xt\‘^mBsV A 

3 ■ 2i07r5 


(B2) 

(B3) 

(B4) 


(B5) 
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Here the dimensionless variables are s = /m^ , rhj = rrii/mB and A = 1 + — 

2 (s + sm\, + rh\-*) and the effective coefficients including the four-quark and gluon dipole 
operators are given by 


= Ct-- 


C3 + - C4 + 206*5 + “^^6 


Cto 

+ — 
4% 


{C,-6C,)A{q^)-CsFP{q^) , (B 6 ) 


= C'9 + h(0,g2) 


- ^h{mb,q‘^) 


■4 11 2 32 ■ 

3 C 1 + 6*2 + y 6*3 - - 6*4 + 526*5 - y 6*6 


[ 4 64 1 

4 

r 16 16 1 

76*3 + -6*4 + 766*5 + y 6*6 

^ 3 

6*3 + y6'5 + y 6*6 


+ 


(y.g 

Att . 

m 

8 -' 

q 


Cl + 4C(<,")) - 3 C 2 ( 2 B (q^) - C (<,=)) - CsCf’ (g=) 


72 


6*1 -l- 776*2^ (1 -l- A„) -|- 26*3 + 206*5 


(B7) 


937 

These include the CKM suppressions and the QCD matching corrections at next-to-leading 
order proportional to A^i = {yubV*g)/{VtbV^g), which corresponds to the small amount of 
CP-violation in the SM. 


Appendix C: Ji coefficients 


In terms of the transversity amplitudes Aq, Ah, A±, and At the Ji coefficients can be 


expressed as [39 


39, 3 


Jl = 


(2 + Pf) 


\Ai\^ + \A\\^ + {L^R) 


+ “7^^® + AfAif*) , (Cl) 


R\2 , 


77 + 2Re {AqA^ ) 


Ji = J [\Al\-‘+ \Al\-‘+ (L ^ R)] , 

J 2 = -jSf [|Aq 1^ + (L ^ R)] , 

J^ = i/3? [|3lit - \Al? + (i ^ R)] . 

J 4 = [Re + (L ^ B)] , 


+ AVsP, 


Js = \/2A 


Re (A^Af) -{L^R)- ^Re (Af A^ + Af A^ 


(C2) 

(C3) 

(C4) 

(C5) 

(C 6 ) 

(C7) 
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J6* = 2A[Re(4^f)-(L^i?)], 
r,=4/3i^Re[A^A*s + {L^R)], 


Jr = V 2 ( 3 i 


Im {AqA^^ ) - (L -)■ i?) + {J^±^*s + J^±^*s) 


Js — ("4o^i ) + (L -)• i?)] , 

J 9 = Pi (^f ^i) + (R —^ -R)] , 


where 


and 




(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 

(C13) 


AA* = Af (g2) Af (g2) + Af (q^) Af (q^) (z, j = 0, ||, ±), (C14) 

in shorthand notation. The transversity amplitndes Af’^ are presented in the appendix for 
A (B) for low (high) region. 

Appendix D: calculation 


The B —)■ K* matrix elements in large recoil limit depend on fonr independent fnnctions 
corresponding to a transversely (a =T) and longitndinally (a =||) polarized K* and at 


next-to-leading order is given by [M| 


Ja — ^aCa + 


Nc rriB 


/ doj /* ^ 

— ^B,±A) / du^K*,a{u)Ta,±{u,u) , (Dl) 

^ Jo 


where S_i_ = 1, Sy = mK*/EK* and the factorization scale qif = ^JmbKqcD- 
The coefficient fnnctions Ca and Ta,± can be written as 

C = 


and 


T - - d I - .] 

/a,± — /a,± H JZ -R:± ■ 


dvr 


The form factor terms at leading order are 




(D2) 


(D3) 


(D4) 
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^(0) _ _ f 


rriB 

2mi) 


y if) 


(D5) 


The coefficients at next-to-leading order can be divided into a factorizable and a non- 
factorizable part as 

C<‘) = cy> + CW) , (D6) 

At NLO the factorizable correction reads 


c5^> = till ^ - i + AmI , 
cj-'* = till^ + 2L + Am\ , 


(D7) 

(D8) 


and the non-factorizable correction for heavy to light transitions are 


fP 


T 


2mhmB 


CaFf) + 2Ci ( Ff ^ + -F 2 


7 ( 9 ) 


( 9 ) 


+ 




i?(7) I r^eff jji(7) I 

GfOii — 02 -r 2 + Os + - - 


8 8 


2mh 


+ 2Fi ( F,'"' + ^F2 


1 ( 9 ) 


( 9 ) 




, (D9) 


where L and AM have given in Ref. Ml- At leading order the hard-spectator scattering 
term (m, u) from weak annihilation diagram is given as 

(m, u) = = 'X 1 + ^) = 0 ’ 

-(0) / .N _ mBU ArriB 


Tn _ {u, oj) = —e. 


^ tubOJ — — ie rrib 


(F3 + 3F4) . 


(DIO) 

(Dll) 


The hard scattering fnnctions at next to leading order contain a factorisable as well as 
non-factorizable part 

P" = ry> + ryfl. (D12) 

Including O (Oi) corrections the factorizable term to the hard scattering functions 7^*2 tire 
given by 

(D13) 
(DM) 
(D15) 


FUu,to) = C‘" 


2 m 


B 


uEk* ’ 


T}_^1 (m, u) = (m, w) = 0 , 


'±, 

-(/) 


TiX {u,u) = 


AmB 


uE 


K* 


and the non-factorizable correction can be compnted by solving the matrix elements of 
fonr-qnark operators and the chromomagnetic dipole operator 

AedCX^ , ms 


rT-(n/) / N 

7]),+ (m,w) = 


+ ^5—M-l (“i ^c) {C 2 + Ca- Cq) 


u + uq^/m\ 2mb 
+ edt± (m, mb) {C 3 + C 4 :-Cq- Amb/mBCXj + (n, 0) Fs] , (D16) 















{ u , w) = 0 , 


(D17) 



= -[e„t|| {u, rric) {C 2 + C 4 - Cg) 

nib ^ 

+ edt\\ {u, mb) {C3 + Ci- Cq) + edt\\ ( m , 0 ) C3] , 


(D18) 


q-(nf) 

'l- 


{u,uj) 


msUJ 


‘^msOJ — — ie 




u + uq'^/m 


B 


6m B 

mb 


h [um% + uq^, me) {C 2 + C 4 + Cq) 


h (um% + uq^, j (C 3 + C 4 + Cq) 
h {um\ + uq^, O) {C 3 + 8 C 4 + 3(5e) — — (773 — C 3 — 15(5g) 


. (D19) 


The ta {u, mg) functions are given by 


^ , q^ 


t± {u,mq) = —^h ("^< 7 ) + - 2^2 (^0 {um\ + uq^,mq) - Bq (g^mg)) , 


t|l {u,mq) = - 


(mq) + (^0 {uml + uq^, mq) - Bq (g^ mq)) , (D20) 




where Bq and Ji are 

Bq [q'^,mq) = —2y^4mq/g2 — 1 arctan 


^ Amy — 1 


Ji imq) = 1 + 


2mA 


f,(rr^2 (^+) + (^+) “ 

a yiiiB q ) 


and 


± T - 




4 + uq"^ 


1 /I 


Li (x) = In - - - In (1 - x) - ^ + Li 2 

X 0 


2 

X 


4 q^ 


X — 1 


(D21) 

(D22) 

(D23) 

(D24) 


Appendix E: Functions involved in Isospin asymmetry parameter 


The function K-) receives an annihilation contribution to leading order and the function 
Ky appear at subleading order of the Ah/mB expansion. The meson photon transition form 
factors and their role in the annihilation contribution to B meson decays has been given by 
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the function with superscript (a). Here only transverse polarization contributes and there 
are no contributions from C 5 and Cq. The function Kf\ contain the decay amplitude from 
the diagram of hard spectator interactions involving the gluonic penguin operator and 
the contributions of the hard spectator interactions diagrams involving the operator (9i_6 
have been included in Kf\. 


The functions K\ and K ^2 dehned via 6 “ is given as 


Kt (,=) = A't'”' (,=)+Ay'"" (y) + (y). 


(El) 


with 


(y) = - (a [pu) + (») . 


1 . 0 ;^* (m) 


Fj_ (s) = - / du 


3 Jq u T us 


1 


"1 


(t>K* (u) 


Ay<‘> (y) = C‘" M (s ), Ai (S) = (S) + 1 / du- 

ms J^c 47r 3 Jq [u + usy 

-4pK.WFv(uml + uq^) (E2) 


^ ^ Nc 471 3 Jo u + us 


and 


with 


{q^) = (g2) + [q^) + [q^) , 

{(f) = + C 2 + ^^4 (h/i) + — 

At"" (fl = 0 ( 41 ) , 

J^±(c) / 2\ Cp O^s {f7h) 1 /" 7 ( (a) / \ 5'_L ^ ('*^) I 77 /- 2 I 2\ 

^2 W ) = -ITT-T—o (“)- 1 - + Mg ) . 


(E3) 


A^c 47r 2 JO 


(E4) 


For parallel case in Eq. (61), 

Kf (q^) = (q^) + (q^) , 


(E5) 


with 


a!<“> (,/) 

a 1"" (y) 

Af"' (y) 


a'A”' (y), 

Cf Qs [l^h) ^ f {um\ + 77 ^) . 




'0 


n + MS' 


iVe 47r 


(E6) 


30 















The vector form factor Fy (s) is given by 

3 

Fy (s) = -{h (s, me, fit) {C2 i/J-h) + C4 (/i^) + Cg M) 

+ h (s, mb, fih) {C3 (fih) + C4 (fih) + Cq {fik)) 

+ h {s, 0, fih) (C's (/^h) + 3^4 if^h) + SCe (/^h)) 

- ^ (Cs M - ft iPk) - 15ft ip,,))} . (E7) 
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